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Abstract 

When the entities undergoing a chemical reaction are not available simultaneously, the classical 
rate equation of a reaction or, alternatively for the evolution of a population, should be extended 
by including non-Mar kovian memory effects. We consider the two cases of an external feedback, 
realized by fixed functions and an internal feedback originated in a self-organized manner by the 
relevant concentration itself. Whereas in the first case the fixed points are not changed, although 
the dynamical process is altered, the second case offers a complete new behaviour, characterized 
by the existence of a time persistent solution. Due to the feedback the reaction may lead to a finite 
concentration in the stationary limit even in case of a single-species pair annihilation A + A ^ 
process. We argue that the different cases are similar to a coupling of additive or multiplicative 
noises in stochastic processes. 
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I. INTRODUCTION 



Chemical reactions involving different species are described by classical rate equations in 
which the time evolution for the forward and reverse reactions are balanced by the product 
of the concentration of the reacting entities, see for a recent review The approach is 
based on the assumption that the reactants are available simultaneously. In the one species- 
annihilation process a particle A is annihilated upon encounter according to the reaction 
scheme A + A ^ 0. The reaction is immediate and is realized with a certain rate u. During 
the related one-species coalescence process A + A ^ A the reactants fuse together realized 
with another rate /i. The time evolution of the global concentration c{t) obeys the mean 
field equation 

^-^ = -{2u-^^y (1) 

whereas u and fi are the coalescence and annihilation rates. In case the reaction takes place 
only after a sufficient accumulation of the reactants, the reaction process is changed crucially. 
We demonstrate that the long-time behaviour is typically dominated by such delay effects. 
The system is able to reach a stationary state instead of exhibiting an algebraic decay in 
time as it follows using Eq. (Q). Let us illustrate the situation in mind in more detail. 
Generally the time evolution of the concentration c{t) is characterized by gain and loss 
terms. Moreover, the concentration could also depend on the history of the sample to which 
it belongs. In the same sense the changing rate of the concentration should be influenced 
by the changing rate in the past. Thus the evolution for c{t) has to be supplemented by a 
memory term indicating the non-Markovian behaviour. Such a term models, for instance the 
way on which a seed concentration had been accumulated by a delayed transport mechanism 
originated by the environment of the reactants. With other words, the changing rate of a 
certain quantity at time t is also determined by the accumulation rate at a former time t' < t. 
In between, i.e. within the interval t = t — t', the reactants are enriched while changing the 
concentration at t'. Regardless that process and further fluctuations the available amount 
of concentration at time t is governed by instantaneous gain and loss of concentration as 
well as on the changing rate at former times t'. Consequently the evolution Eq. (P) should 
be modified according to 

dtc{t) = rc{t) - uc\t) - [ K{t,t';c)dt'c{t')dt' . (2) 
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The last term represents a permanent re-evaluation of the reactions. Hence the aim of 
the present paper is to study models which follow an evolution equation as given by Eq. ^ 
where the memory kernel K will be specified below. As a new ingredient we assume that the 
memory term may depend on the concentration c(t) and its derivative as already indicated in 
Eq. Whereas the first realization consists of a kernel, given by a deterministic function, 
the second case is characterized by a kernel, determined by the concentration itself. The 
first case is denoted as an external feedback with a distributed time-delay, the second one is 
signified as an internal delay, since the time scale of the kernel is defined in a self-organized 
manner by the time scale of the concentration itself. 

Our model can be grouped into the increasing interest of incorporating delay and feedback 
mechanism in a large variety of systems. Obviously the crucial factor governing the dynamics 
of systems comprising many "units" consists of interaction and competition. Such features 
are believed to underlie the complex dynamics observed in disciplines as diverse as economics 
P], biology and weather j^, politics jj], medical care j^, and ecology 0]. As discussed 
recently and also stressed by the present paper another characteristic trait of physical 0, 
as well as biological systems is played by the mentioned time-delayed couplings. Such 
memory effects, considered in a large variety of systems, could be a further unifying feature 
of complex systems. It is well-known that evolution equations with memory kernels as it is 
presented by Eq. (|21), can be derived following the well established projector formalism due 
to [l^, see also [H]. That approach had been applied successfully for the density- density 
correlation function in studying the processes in undercooled liquids 3, ll^- Recently 
a Fokker-Planck equation with a non-linear memory term was used to discuss anomalous 
diffusion in disordered systems ilM- The results could be confirmed by numerical simulations 



including diffusion on fractals 



ISj. see als o (igI . Ii7 | . Moreover, it had been demonstrated 
that mobile particles remain localized [l8| due to the feedback-coupling. Notice that a 
formal solution of whole class of non-Markovian Fokker-Planck equations can be expressed 
through the solution of the Markovian equation with the same Fokker-Planck operator . 
The non-Gaussian fluctuations of the asset price can be also traced back to memory effects 
I20I 1. An additional cumulative feedback coupling within the Lotka-Volterra model, which 



may stem 
behaviour 



Tom mutations of the species or a climate change, leads to a complete different 



21[ compared to the conventional model. If the Ginzburg-Landau model for the 



time evolution of an order parameter is supplemented by a competing memory term, the 
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asymptotic behaviour and the phase diagram is completely dominated by such a term 22 1. 



Different to our self organized approach with K{t, t'; c) there is a broad class of models 



with delay integrals without a dependence on the variable c{t) [23|, |2J, |25[ , for a survey 
and applications in biology see |26|. The spreading of an agent in a medium with long-time 
memory, which can model epidemics, is studied in ^7}. Time-delayed feedback control is an 
efficient method for stabilizing unstable periodic orbits of chaotic systems 28|. It may also 
induce various patterns including travelling rolls, spirals and other patterns 29]. A further 



3- 



approach concerns the study of bistable time-delayed feedback systems driven by noise 
Based on a former paper the distribution function of the first-passage-time is analysed 
where both, the delay and the noise term offers a significant influence on the behaviour. A 
global feedback is studied recently also in a bistable system . The purpose of that paper 
is a discussion of the domain-size control by a feedback. Even in an open quantum system 
non-Markovian dynamics is characterized by a time-non-locality in the equation of motion 
for the reduced density operator j^^. 

In view of the large variety of systems with feedback couplings it seems to be worth to study 
simple models, which still conserve the crucial dynamical features of evolution models as 
non-linearities and moreover as a new ingredient, delayed feedback-coupling. In the present 
paper we demonstrate the crucial influence of the non-Markovian memory term on chemical 
reactions, where the retardation effects are comprised into a memory kernel K{t). That 
memory kernel yields an additional competitive term to the instantaneous non-linear terms. 

II. EXTERNAL FEEDBACK 
A. Discrete time-delay 

In this section we specify the model, Eq. Q, by fixing the memory kernel K in terms 
of deterministic functions. This procedure is well-known in the framework of population 



dynamics (an introduction is given in j3^). Different to these approach our rate equation 
(0) offers a feedback-coupling to the decay rate of the population or the density of a chemical 
species dtc(t), respectively. By choosing the kernel K(t,t'; r) = fi6(t — t' — r), and substi- 
tuting K in the dimensionless form of Eq. Q, one gets a discrete time-delayed differential 
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equation 

^-^ ^ c(t) - cHt) - . (3) 

Here the quantity r is the delay time which is assumed to be short ranged t <^t. If r = 0, it 
results in a modification of the evolution equation for the logistic growth, which is solvable 
(see below). The limiting case r ^ can be obtained by expanding the delayed term of Eq. 
0. This leads to an equation 

(Pc 1 + li dc 1 , , , , 

^- + — cl-c=0 /i^O, 4 

dt^ dt fiT 

which cannot be solved in closed form due to the non-linearity. One gets an analytical 
solution, if one neglects the non-linearity, which corresponds to the exponential growth of 
the particle density in addition of a feedback-coupling. Using standard methods the solution 
c(t) is given by 

c(t) = Ci exp(K+t)+C2 exp(K_ t) with k± = (l + ^± + /i)^ — 4/ir) . (5) 

2 flT \ ) 

The integration constants can be calculated to C\ = cq '^'^'^ and C2 = cq , if one 

suggests c(0) = Co and ^ = cq. Discussing the stability of Eq. (0) one finds out, that for 
fi > the both exponents k±. are positive. In case of < one of them is positive, i. e. there 
are no bounded solution in the limit t 00. Alternatively, one can treat the differential 
equation with delay by using the Laplace transformation. Denoting the transformed function 
by c{z) = C{c(t)} = c(t) exp{—zt)dt the transformed function of the linear equivalent of 
Eq. obeys the relation 



Co [! + (! + /u) exp(-r 2;)] - c(-r) + z exp(-r 2;) / exp(-r^) c(^) d^ 

c(z) = : -, . (6) 

^ ^ z [1 + fi exp{-T z)] - 1 ^ ^ 

As displayed by the last relation, differential equations with delayed terms are distinguished 
from conventional equations by one important feature, that the function c{t) has to be given 
within the whole interval — r < t < and not only at one special point c{t = 0) = Cq. Fixing 
c{t) = Co in that interval — r < t < 0, then the Eq. (jH)) becomes simpler 

ciz) = ^ "\ — . (7) 

l+fl Oxp( — T z) 

This equation has to be transformed back. Since the zeros of the denominator of Eq. (|7j) are 
not available analjd;ically, the function c(t) could not be not calculated exactly. However, 



one can convince that in the long time hmit, ^ 1, the concentration c{t) offers an 
exponential form, already obtained in Eq. (|^. 



B. Distributed time-delay - Exponential kernel 

An usual choice of a kernel in problems dealing with continuously distributed time-delay is 
the exponential kernel K{t, t') = fi exp [—A {t — t')] (A > 0) , which satisfies some properties 
like the boundedness and positivity of the kernel. The parameter A determines the time 
scale of the memory. Substituting this kernel in dimensionless form of Eq. (0) it leads to 
the following equation 

t 

dtc{t) = c(t) - c^{t) - /i y exp [-A {t - t')] dtic{t') dt' c(0) = cq . (8) 



Fixed points of this equations are Cg = and = 1. They are independent of the memory 
parameter fi which is a general result, compare the discussion at the end of this section. To 
gain information about their stability, we investigate the linear equation in analogy to Eq. 
(jS)) firstly. Using Laplace transformation results in 

[z - 1) [z + X) + fiz 



The zeros of the denominator of the latter formula are given by z± = 1/2 [— (A + yU — 1)± vD] 
with D = {\ + fi — 1)"^ + 4 \. It is obvious that D > 0, because the assumption A > 0. 
Further one observes that > |A + /i — 1| and so z+ is positive and z_ is negative. Because 
the structure of the solution is c(t) = A exp{z^ t) + B exp(z_ t) , we found only unstable or 
unbounded solutions (for t — > oo), respectively. That means that the exponential growth, 
which is observed in case of = 0, could not be restricted through the additional feedback 
term. Alternatively, Eq. (jH)) can be analysed by a repeated differentiating with respect to 
t and a renewed use of the first order equation to eliminate the integral. The procedure 
results in the second order equation 

^ + [2c + A + /i-l]^ + Ac[c-l] = with c(0) = Co and = cq [1 - Cq] (10) 

dt^ dt dt 

that linear equivalent exhibits the same behaviour found above. The non-linear Eq. 
gives rise to a more complex stability behaviour. To get insight into the stability behaviour 
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of both fixed points we take tlie ansatz c(t) = Cs + ^{t). Substituting tliis ansatz into Eq. (jSl) 
and omitting all terms of order ip^ one gets the equations for Lp[t) 

t 

>p{t) = ± ^(t) -fi y exp [-A (t - t')] <p{t') dt' , (11) 







where the positive sign corresponds to = and the negative one to = 1. The equation 
for the positive sign is already discussed above, and it follows that the trivial solution 
Cs = is an unstable fixed point. To gain the stability behaviour of the second fixed 
point we follow the same steps as above. The Laplace transform of (f{t) is denoted as 
if{z) = C{f{t)}{z) = (f{t)exp{-zt)dt. We get 



In this case the zeros of the denominator are given by 1/2 [— (1 + A + /x) ± vD], where the 
discriminate is D = (1 + A + fi)"^ — 4 A. The solution of (p{t) (for D ^ 0) can be written as 
ip{t) = (fo [A exp(z+t) + B exp{z^ t)] with the coefficients 

^^Zl±A±^^±^ and g,l-(A + AO + v^ ,j3) 



The last formula implies that for D > the solution is real and hence physically relevant. 
The limiting case D = delivers critical values either for /i or for A. Choosing [i = /i(A) the 
critical values for /i are /i^(A) = -(A + 1) ± 2 ^A. If yU is situated in the intervals ^ > ^\ 
or fi < fit, D > is fulfilled and the resulting solution is real. The solution for (f is only 
stable, if /i > — (1 + A). All together the parameter has to fulfil fi > fi\ to get stable, 
physically relevant solutions. In the limiting case, fi = fJi\, one specifies 



(/p(t) = (/?o exp (^-v^t) l + {\/\-l)t , (14) 



which is a bounded solution for t — oo. In case oi \ < 1, (p{t) changes its sign at t* = 
(1 — \/A)^^. Together we can determine the stability domain for the fixed point = 1 to 
the area in the (A, /i)-plane, where > /i+(A) is fulfilled. Eq. (fTIHl is an evolution equation, 
which is well-known as equation for a damped oscillator. Interpreting the equation as an 
equation of motion for a particle in a potential U (c) = — c^/2 + then the factor in front 
of the first derivative has the meaning of the damping constant. Due to the non-linearity 
of the underlying evolution equation the damping parameter 7(c) = 2c{t) + A + — 1 is 



driven by the time dependent concentration. If one regognizes that only a positive damping 
parameter 7(c) > is reasonable, it is obvious that the stability criteria is additionally 
depended on the initial value Cq. That is indeed a feature of systems including feedback or 
memory effects. With the following constraints 

fi + X + l-2y/X >0 
2co + /x + A-l > (15) 

one gets the stability domain depicted in FigU where the (/i, A, Cq)- plane is displayed. The 
area, limited by the curved area and the plane area at the forefront, is the region where the 
solution is supposed to be stable. 



C. Distributed time-delay - a periodic kernel 

Since the dynamical behaviour is affected by the sign of the memory parameter /i, we 
consider now a periodic kernel such as K{t, t') = fj, cos [A {t — t')], which results in 

t 

dtc{t) = c(t) -c\t)~fi J cos [A (t - t')] dt>c{t') dt' . (16) 



Again the two stationary values are determined to = and Cg = 1. The Laplace transform 
of the corresponding linear evolution equation is calculated to 

The denominator is a cubic polynom and the zeros could be calculated. At least one of the 
three solutions has a positive real part in case of A 7^ and so there are only unbounded 
solutions for the linear equivalent of Eq. (fTBj) for long times. Like in the case of the expo- 
nential kernel Eq. ()16|) can be written as a conventional differential equation which is in the 
present case third order 



d^c , , d'^c 

_+[2c + ^-l]- + 



dC , n 



| + A'^c(c-1) = (18) 



with the initial value c(0) = cq, the initial slope = Cq (1 — Cq) and the initial curvature 



Co (1 — Co) (1 — yU — 2 Co). The invar iance of the kernel against the inversion of A ^ — A 
is maintained in Eq. (fTH|). too. The third order differential equation is again equivalent to a 



system of three first order equations. The linear form leads to the same results as discussed 
after Eq. (fTTj) . In FigI21we show a typical example for the time evolution of the concentration. 
Obviously, an increasing value of the parameter A leads to an overdamped behaviour. 
The special case of a constant memory kernel K{t, t') = fi can be treated exactly. Setting 
A = in Eq. (fTB|) it results the following equation 

dc 

— = [l-fi]c{t)~c\t)+fiCo . (19) 

The decay of the concentration decay depends directly on the initial value cq making the 
memory apparently. Further the growth rate is modified by replacing 1 — > 1 — /i. The 
general solution for arbitrary fi is deducible with the result 



c(t) 



1 



D 



tanh 



Dt + c 



with D = a/(1 — + 4/i Co and c = - In 



D + 2co + fi-l 



(20) 



D - 2co - /i + 1 

The last solution is only reasonable within the interval < cq < 1. The special case fi = 1 
yields 



c(t) 



Co tanh 



Vc^t+ - In 



1 



(21) 



1 - Vco^ 

The characteristic width of the curve To = 2D^^ is determined by the memory parameter n 
and the initial value Co. For t oc the solution, given by Eqs. (HHIEII), lead to fixed points 

C,(/i) = ^[l-/i±D] 

different from c^ = or c^ = 1, respectively. The negative branch is omitted, because c^ 
has to be positive. Contrary to the previous cases the stationary solution depends on the 
memory parameter /i. Concerning the non-trivial case, where A 7^ 0, the instability of 
Cs = is maintained for all values of /i, whereas c^ = 1 is stable in the half-plane fi > 0. In 
this area one gets a solution starting with an exponential decay (growth) which is followed 
by a damped oscillational motion toward the stationary solution c^ = 1, compare Fig|21 



The analysis can be generalized to an arbitrary kernel K{t). Considering the Laplace 
transformed of Eq. Q, one concludes straightforwardly, that the fixed points remain un- 
changed, whenever the transformed memory kernel satisfies the relation 

\imzK{z) = . (22) 

2:— >o 
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For a constant kernel that equation is not fulfilled and one finds the modified fixed point 
given above. Thus, an external non-trivial feedback-kernel behaves like an additive noise 
in stochastic processes. In the next section we demonstrate how that picture is altered by 
assuming an internal feedback. In that case the kernel is determined by concentration itself. 
As consequence the last relation is likewise violated and the stationary points are directly 
controlled by the memory. Such a situation is the feature of stochastic processes whenever 
a multiplicative noise is included. 

III. INTERNAL FEEDBACK 

In this section let us generalize the previous approach by including self-organized delay 
effects where, as already discussed in the introduction and in accordance with the conclusion 
of the last section, the kernel depends on the concentration c(t) itself. That means, the time 
scale of K{t) in Eq. is determined by the time scale of c{t). Following this line that the 
delay effects are dominated by the concentration c{t) we conclude 



The kernel is defined by the concentration in the intermediate interval t — t' which is coupled 
to the changing rate at the previous time t' as indicated in Eq. Q. The memory kernel 
characterizes the way on which a seed concentration had been accumulated. Within the time 
interval t = t — t', the concentration is further enriched, or with other words, at time t', 
only an incomplete reaction is realized. During the time interval t — t' the residual particles 
are moved to the reaction zone. As the simplest realization of Eq. fj23|) we chose a linear 
dependence on the concentration. In that case the memory term is a competitive one to the 
conventional quadratic term. Summarizing all contributions we will analyse the evolution 
equation 



The first term characterizes the spontaneous creation of particles (with rate r > 0) in 
according to the reaction ^ A, whereas the second one describes a two particle reaction 
discussed in the introduction, compare Eq. (^. Both terms leads to a stable stationary 
solution when they are balanced, manifested by r > 0, m > 0. For generality one could 
also discuss the case that the parameter r and u can vary and change their signs. Here, 




(23) 




(24) 
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we restrict us to positive quantities r and u. The last term models the feedback with the 
coupling parameter k which mimics the influence of the environment. The following results 
are strongly influenced by the sign of the memory strength k. In case of k < and m > 
both non-linear terms are competitive ones. Let us remark that the memory kernel gives rise 
to a coupling of the time scales. In the vicinity of the upper limit of the integral t' ~ t the 
memory term reads c(0)9tc(t), i.e. a momentary change at the observation time t is coupled 
to the value at the initial time t = 0. Therefore the very past is related to the instantaneous 
value of c{t). In the opposite case, at the lower limit t' ^ 0, the change of the concentration 
near to the initial value df c{t' = 0) is directly coupled to the instantaneous value c{t). In 
such a manner the memory term represents a weighted coupling of the behaviour at the 
initial time and the observation time. This coupling leads to another long-time behaviour. 
Notice that the generic behaviour, discussed below, is not changed by assuming other terms 
with different power laws in Eq. ()24|) . 

A. Stationary solutions 

In this section we find the stationary solutions of Eq. ()24|) and discuss their stability. 
An important case is realized by r > 0, m > and arbitrary k. Introducing dimensionless 
quantities c r/uc, t ^ t/r we end up with an equation 

dtc{t) = c{t) - c\t) -11 f c{t- f) dt>c{t')dt' (25) 

with = k/u. The solution of the evolution Eq. is simple when the memory kernel is 
zero, i.e. /i = 0. It results 

^^^) = l + co(e*-l) ^'^^ co = c(t = 0) . (26) 
In the long time limit it results a non-trivial stationary and stable solution Cs{fi = 0) = 1 
and an unstable trivial solution Cg = 0. The inclusion of a memory term will change that 
behaviour drastically. To that aim let us analyse how the concentration is controlled by the 
feedback-coupling strength fi. The formal solution of Eq. is obtained by using Laplace 
transformation. We get 

_ co{l + fic{z)) -A{z) 
"^^^ ~ z{l + ^Iciz)) + l 

with A{z) = C{c\t)){z) . (27) 
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Remark tat we discuss only the non-trivial memory controlled solution by assuming cq ^ 
and Co 7^ 1. Following the line of the previous section the long-time behaviour is obtained 
by making the ansatz c{t) = Cs + if{t) or after Laplace transformation 

c{z) = ^ + ^{z) (28) 

where the function ip{z) remains regular for z ^ 0. This function will be discussed in the 
next subsection in combination with the linear stability. The quantity Cs(/i) represents, like 
before, the stationary solution in the limit t oo. Apart from the trivial solution Cg = 
there exists a non-trivial, memory controlled solution 

1 + /iCo 

Csifi) = ■ (29) 

1 ~r 

Notice that this result is in accordancce with the general conclusion expressed by Eq. (|22j). 
The stationary solution, Eq. (jSHI), depends on both, the initial value cq and the memory 
strength /x. If one applies the Laplace transformation in Eq. Q and making linear stability 
analysis, see the subsequent subsection one ends up only with the memory independent fixed 
points Cg = 1 (stable) and Cg = (unstable), respectively. Here, we have demonstrated that 
an internal feedback leads to a modified non-trivial stationary fixed point. Because c{t) is a 
concentration the stationary solution is only accessible in case Cg > 0. The non-Markovian 
behaviour is apparent by the dependence of Cg{fi) on the initial concentration cq and on the 
memory parameter /i. 



B. Phase diagram 

In this subsection the phase diagram is obtained employing a linear stability analysis. 
Notice that it can not be performed in terms of y^{z) defined in Eq. p8|l . Instead of that we 
have to insert c(t) = Cg{fi) + ip{t) in Eq. (j^ . As the result one finds ip oc exp(— At), where 
the stability exponent A reads 

A = -1 if c, = 

^^c^l + 2c^ if . (30) 

X ~|~ ^ 

The phase diagram is defined in the /x—co— plane under exclusion of the trivial case cq = 1. In 
the stable region both conditions, > and A > 0, has to be fulfilled simultaneously. Notice 
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that there exits regions where the stationary solution is available Cg > but not stable A < 0. 
The different cases are summarized in Fig|21 The non-trivial stationary solution is stable for 
/i > — 1 in the interval < Cq < 1 and for cq > 1 in case of > /xi (cq) = — 1 + a/1 — Cq ^ 
There exists a second domain where the non-trivial stationary solution should be stable. 
This area is limited by the lines fj, = —1 and /i2(co) = — 1 — a/1 — c^^ (dashed lines in 
FigO)). However the system can not achieve this region, because the initial value cq > 1 
leads to Cs(yu) > Cq in contradiction to the monotone decrease of the function c(t) as indicated 
by Eq. ()25|). We omit the discussion of the other cases with different signs of the parameter 
r and u in Eq. ()24|) . The inclusion of an internal memory leads to a non-trivial stationary 
point which is controlled strongly by the memory parameter fi. Because the feedback is 
characterized by the concentration, relation Eq. (j^^ is violated leading to a new fixed point 
in comparison to the non-Markovian situation. 



IV. CONCLUSIONS 



In the present paper we have generalized the conventional rate equations of the reaction 
limited case by including memory effects. Once, the memory can be originated by external 
constraints which hints the particle to encounter simultaneously. As a consequence the rate 
of the concentration at the present time may be coupled to the rate at a previous time, 
where between both processes an additional delay-time r appears. Although the stationary 
point are unchanged by such kind of external delay, the dynamical behaviour to reach that 
fixed points is altered. This situation reminds to an additive noise coupling in stochastic 
processes. In case of an internal feedback the stationary solution is changed and depends on 
the memory coupling and the initial value. The reason for that consists of the self-organized 
manner inherent in the memory. The time scale of the relevant variable, the concentration, 
determines the time scale of the feedback. Both realizations are characterized by an explicit 
coupling of the rates at the observation time t and a previous time < t' < t or with other 
words, the instantaneous changing rate of the concentration is not only determined by time- 
local gain and loss terms, but additionally by the changing rate in the past. Such a situation 
may realized whenever the reacting entities are not available simultaneously. Thus, the 
reaction is determined by the accumulation of the species at previous times. To capture the 
influence of that delay process, a memory term is included into the evolution equation. The 
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non-Markovian part is generally found adopting projection methods of statistical mechanics. 
Especially, the form of the memory is suggested by analyzing glasses or anomalous diffusion 
in disordered media. Although we are aware, that the analytical expression of the feedback 
term is the most controversial point, the delay is in according to the general approach 
applied for chemical reactions. The memory term allows to study a real competitive situation 
between two non-linear terms in our evolution equation, namely the instantaneous loss term 
and the memory driven one. Even that kind of competition leads to the richer behaviour 
discussed in the present paper. Due to the memory the present and the initial times are 
mixed. Therefore, the system offers the possibihty that the reaction process can be reduced 
or or speeded up. As the result wc get a new stationary behaviour which differs significantly 
from the standard one. We believe that memory effects are also a feature of other dynamical 
complex systems. 
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FIG. 1: Stability domain for an exponential kernel in the (/x, X,cq) plane. 




FIG. 2: Time evolution of c{t) with a periodic kernel and A = 0.1, /U = 1 and cq = 0.5. 
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FIG. 3: Stability of the internal memory controlled stationary solution in the cq — //-plane. 
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